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We present a measurement of Z boson pair production in pp collisions at 1.96 TeV with 9.6 fb^^ 
to 9.8 fb^^ of DO data. We examine the final states eeee, eefifi, and fifififi. Based on se- 
lected data, the measured cross section in the mass region M{Z/'y*) > 30 GeV is a-{pp — >■ 



1.05: 
the ZZ 



Zj^* Z/'y*) — 1.26_g'3g (stat)_Q\5 (syst) ± 0.08 (lumi) pb; after correcting for the expected ra- 
tio of (j{pp — >■ Z/y* Z/y*) to (j{pp -^ ZZ), we derive a cross section for pp — >■ ZZ production of 
-0 30 (st8't)_o 12 (syst) ± 0.06 (lumi) pb. This result is combined with a previous result from 
t^H.^vv channel resulting in a combined pp — > ZZ cross section measurement of 
1.32tQ'25 (stat) ± 0.12 (syst) ± 0.04 (lumi) pb. These measurements are consistent with the stan- 
dard model expectation of 1.43 ± 0.10 pb. We extend this analysis to search for the standard model 
(SM) Higgs boson between 115 and 200 GeV. At a Higgs boson mass of 125 GeV, we expect to set 
a limit of 43 times the SM expectation at 95% C.L., and set a limit of 42 times the SM expectation 
at 95% C.L. 



PACS numbers: 12.15.Ji,13.85.Qk,14.70.Hp,14.80.Bn 



I. INTRODUCTION 

We present a measurement of the cross section a{pp — > 
Z/7* Z/y*) at ^/s — 1.96 TeV, using events where each 
Z/7* results in two charged leptons. Because the branch- 
ing fraction of the Z boson to charged leptons is smaller 
than that to quarks or neutrinos, this process is relatively 
rare, but has the advantage of being an extremely pure 
final state. The largest fraction of the background results 
from events in which one or more jet has been misiden- 
tified as a lepton, since few other processes in the stan- 
dard model (SM) produce four isolated leptons. We also 
unfold our measurement to determine the crijap — > ZZ) 
cross section. 

After measuring the i-channel Z/y* Z/y* cross sec- 
tion, we reinterpret the analysis as a search for the Higgs 
boson in the four lepton final state, predicted in the SM 
as a result of electroweak symmetry breaking. Both the 
ATLAS and CMS experiments at the CERN LHC pp col- 
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lider have observed a four lepton resonance at a mass of 
~^125 GeV [1, 2] which, when combined with other decay 
channels, is consistent with the SM Higgs boson. 

Z boson pair production was studied at the CERN 
LEP2 colhder by the ALEPH [3], DELPHI [4], L3 [5], 
and OPAL [6] collaborations in multiple final states, in- 
cluding e+e^ — )■ i^i~(. +£ ~, where £ represents an elec- 
tron or a muon. The LEP experiments also set limits on 
anomalous ZZZ and ZZy couplings [7]. 

The Fermilab Tevatron experiments have also searched 
for and measured the pair production of Z bosons. The 
DO collaboration's analysis of ZZ -^ £^£^£ +i? ^ produc- 
tion with 1.1 fb^^ oi pp data yielded an upper limit of 
4.4 pb on the ZZ production cross section at 95% C.L. 
Additionally, limits on anomalous ZZZ and ZZj cou- 
plings were determined [8]. The DO collaboration was 
the first to observe ZZ production in pp collisions in the 
£+£-£'+£'- final state with 2.7 ib'^ of data [9]. The DO 
collaboration has also measured the ZZ cross section in 
the £'^£^vv final state, first with 2.2 fb^^ [10] and later 
with 8.6 fb""'^ of integrated luminosity, yielding a final 
measurement of 1.64±0.44(stat)+g}5 (syst) pb [11]. The 
CDF collaboration has analyzed data from 1.9 fb^^ of 
integrated luminosity to study ZZ production, measur- 
ing, when combining £'^£~£'^£ ~ and £'^£'vv channels, 
a cross section of a{ZZ) = lAt_QQ (stat + syst) pb [12]. 
The ATLAS collaboration has observed pp — )■ ZZ pro- 
duction in the four charged lepton final state in 1.0 fb^^ 
of data at Vs = 7 TeV [13]. The CMS collaboration has 
measured c!{pp -^ ZZ) in 5.0 fb~^ of data at -y/i = 7 TeV 



[14], and has observed the rare decay Z — > i'^i~i ^i ~ 
with a branching fraction in agreement with the SM pre- 
diction. 

This article is an update of the DO cohaboration's prior 
ZZ to four charged lepton analysis that measured a cross 
section of cr(pp -^ ZZ) = 1.26to:3^ (stat) ± 0.11 (syst) ± 
0.08 (lumi) pb using 6.4 fb~^ of integrated luminosity 
[16]. The result presented here uses 9.6 fb^^ to 9.8 fb^^ of 
integrated luminosity, and expands electron acceptance 
in the eeee final state. 



II. DETECTOR 

The DO detector is described in detail elsewhere [17- 
20]. The main components are the central tracking 
system, the calorimeter system, and the muon detec- 
tors. The central-tracking system is located within a 
2 T solenoidal field and consists of two different track- 
ers. Located closest to the interaction point is the silicon 
microstrip tracker (SMT) and surrounding that is the 
central fiber tracker (CFT). The SMT is an assembly of 
barrel silicon detectors in the central region, along with 
large-diameter disks in the forward regions for tracking at 
high pseudorapidity (77) [21]. The CFT consists of eight 
concentric coaxial barrels each carrying two doublet lay- 
ers of scintillating fibers. The liquid-argon calorimeter 
system is housed in three cryostats. The central calorime- 
ter (CC) covers up to \ri\ = 1, and two end calorime- 
ters (EC) are located in the forward regions, extending 
coverage to \ri\ = 4. In the intercryostat region (ICR) 
between the CC and EC cryostats, there is a scintillat- 
ing intercryostat detector (ICD) between 1.1 < \ri\ < 1.4 
that recovers some energy from particles passing through 
the ICR. Closest to the collisions are the electromagnetic 
(EM) regions of the calorimeter followed by hadronic lay- 
ers of fine and coarse segmentation. 

A muon detection system [22] is located beyond the 
calorimeters and consists of a layer of tracking detectors 
and scintillation trigger counters before 1.8 T toroid mag- 
nets, followed by two similar layers after the toroids. 

There is a three-level trigger system consisting of a 
collection of specialized hardware elements, microproces- 
sors, and decision-making algorithms to selectively record 
the events of most interest. 



III. MONTE CARLO 

We use the pythia [23] Monte Carlo (MC) program 
to determine the Z/7* Z/j* -^ £+£^£'+£'^ signal accep- 
tance and to simulate the migration background. The 
signal is defined to consist of Z/j* Z/'y* pairs where each 
Z/j* boson has a mass greater than 30 GeV. The mi- 
gration background consists of Z/'y* Z/y* events where 
at least one of the two Z/y* bosons has an invariant 
mass of less than 30 GeV; it enters the signal sample 
either due to mismeasurement or by mis-assigning the 



lepton pairs in the eeee and ij-ij-ij-ij- channels. We in- 
clude Z/y* Z/j* — > £+£^T^T^ events where the taus 
decay into electrons or muons as appropriate to match 
the final four-lepton signature in the signal acceptance. 
Contributions from ZZ — >■ t^t~t^t~ with subsequent 
decays into muons and electrons are also examined, but 
found to be negligible. The ZZ transverse momentum 
(pt) spectrum is also estimated using SHERPA MC [24], 
and the difference between the pr spectra from pythia 
and SHERPA is used as a systematic. The dominant tree- 
level diagrams for pp — )• Z/y* Z/y* -^ £+£^/+£'~ are 
shown in Fig. 1. The singly resonant Z boson diagram 
contributes at low mass, and we expect a negligible con- 
tribution to the signal yields from this diagram in our 
analysis. 




FIG. 1: Feynman diagrams for (a) the t-channel tree-level 
process qq — > ZZ — > £'^£^£ ^£ ^ and (b) the singly resonant 
process. 

To estimate the migration background, we generate 
Z/y* pairs where at least one of the bosons has a mass 
between 5 and 30 GcV, and estimate the cross section of 
these events using next-to-leading-order (NLO) MC from 
MCFM [25] with the CTEQ61M PDF set [26]. 

The ti background is estimated using ALPGEN [30] with 
a top quark mass of 172 GeV and is normalized to an 
approximate NNLO cross section calculation [31]. 

Other backgrounds where photons or jets must be 
misidentified for the event to enter our sample, such as 
events containing a Z plus jets, are estimated from data 
as described in Sec. VI. 

For the Higgs boson search, we generate SM Higgs bo- 
son events with masses between 115 and 200 GeV in 5 
GeV increments. We simulate the gluon fusion {gg -^ H) 
and ZH associated production [qq -^ ZH) processes 
using PYTHIA. The expected gg ^ H cross section is 
corrected to next-to-NLO (NNLO) with next-to-next-to- 
leading- log resummation of soft gluons [27]. The associ- 
ated ZH production cross section is corrected to NNLO 
[28]. The expected branching fractions for the Higgs bo- 
son decay are determined using hdecay [29]. 

All of the MC samples are passed through a GEANT [32] 
simulation of the DO detector. To account for detector 
noise and additional pp interactions, data from random 
beam crossings are overlaid onto all MC events to match 
the instantaneous luminosity distribution of the selected 



data. The same algorithms used to reconstruct real data 
events are run on these simulated events. 



IV. OBJECT IDENTIFICATION 

All muon candidates are reconstructed either as a 
muon track reconstructed from hits in both the wire 
chambers and scintillators in at least one layer of the 
muon system, or as a narrow energy deposit in the 
calorimeter system consistent with that expected from 
a muon passing through the calorimeter that is not asso- 
ciated with tracks in the muon system. Each muon can- 
didate must be matched to a track in the central tracker 
with a pt > 15 GeV, and the track pT is taken as the 
Pt of the muon, pi^. This track must have an impact pa- 
rameter consistent with the muon coming from the inter- 
action point. We consider two muon isolation variables: 
, the scalar sum of the track pT within a cone of 



TTitrkconc 

H/rp 



AR < 0.5 [33] about the muon track; and £'^^'°, the sum 
of the calorimeter energy in an annulus 0.1 < AR < 0.4 
centered on the muon track. If the muon is reconstructed 
in the muon system, then we impose the requirement that 



trkcono /^M < 0.25 and EY°/Pt < 0-4- Otherwise, each 



^trkcono/p; 



variable divided by pi^ must be less than 0.1. 

Different selection requirements apply for electrons 
identified in the CC (|r/d| < 1.1), EC (1.5 < \r]d\ < 3.2), 
and ICR (1.1 < \r]d\ < 1.5), where rjd is the pseudorapid- 
ity calculated with respect to the center of the detector. 
In the CC and EC, electrons must have at least 90% of 
their energy found in the EM calorimeter, have p^ > 15 
GeV, and pass a calorimeter isolation requirement. The 
Pt estimate for the CC and EC electrons is based on the 
energy deposited in the calorimeter. For electrons in the 
CC, the sum of transverse momenta of the charged cen- 
tral tracks in an annulus of 0.05 < AR < 0.4 about the 
electron, 74, must be less than 4.0 GeV. There must ei- 
ther be a track in the central tracker associated with the 
calorimeter cluster, or hits in the central tracker consis- 
tent with a track along the extrapolation of the calorime- 
ter cluster to the interaction point. Finally, the electron 
must pass a neural net (NN) discriminant trained to sep- 
arate electrons from jets in the CC using seven shower 
shape and isolation variables as input. 

In the EC only, we require that the track isolation 74 
be less than (7.0-2.5x |r/d|) GeV or 0.01 GeV, whichever 
is larger. The electron must pass a NN discriminant 
trained to separate electrons from jets in the EC using 
three shower shape and isolation variables as input and 
an additional chi-square-based shower shape requirement 
designed to distinguish electrons from jets. 

Within the ICR, there is incomplete EM calorimeter 
coverage, so the electron must pass a minimum EM + 
ICD energy fraction requirement that varies with \rid\. 
The candidate must be matched to a central track with 
Pt > 15 GeV and have a pt > 10 GeV measured in the 
calorimeter. Additionally, the ICR electron must sat- 
isfy two multivariate discriminants designed to reject jet 



background. Due to the limited energy resolution in the 
ICR, we use the pr of the track associated with the ICR 
electron to estimate the ICR electron energy. 

Jets are used in the estimation of the instrumental 
background, as discussed in Sec. VI. In this analysis, we 
use jets reconstructed from energy deposits in the CC, 
EC, and ICD detectors using the Run II midpoint cone 
algorithm [34] with a cone size of AR — 0.5. The jets 
must have px > 15 GeV and \i]d\ < 3.2. We apply the 
standard jet energy scale (JES) corrections [35] to jets in 
both data and MC. 

The missing transverse energy, _^j,, is calculated using 
a vector sum of the transverse components of calorime- 
ter energy depositions, with appropriate JES corrections 
[35]. In the eefifi and fJ.fJ,fJ.fJ- final states, the $!rp is cor- 
rected for identified muons. 



V. EVENT SELECTION 

To maximize the acceptance, we consider all events 
that pass the event selection requirements below without 
requiring a specific trigger. The majority of our accep- 
tance comes from events collected by single lepton and 
di-lcpton triggers. As there are four high-p^ leptons in 
this final state, we estimate that the trigger efficiency for 
the signal is greater than 99.5% in all channels. 



A. eeee final state 

All electron candidates have to satisfy the requirements 
in Sec. IV. We require at least four electron candidates. 
If there are four CC/EC electron candidates, no ICR elec- 
tron candidates are considered, and if there are more than 
four CC/EC electron candidates, the highest-pT candi- 
dates are used. At least two of the electrons must be in 
the CC, and if an event has more than one ICR electron, 
only the leading ICR electron is considered as a lepton 
candidate. All possible pairings of the selected electrons 
are considered with no charge requirement imposed, and 
we require that one of the pairings has di-electron mass 
Mee > 30 GeV for both di-electrons. Additionally, there 
must be AT? > 0.5 between any ICR electron and any CC 
and EC electrons, or the ICR electron is not considered. 
Because the instrumental background contamination is 
expected to vary significantly depending on the number 
of central electrons, the eeee channel is then divided into 
four sub-channels that depend on the number of elec- 
trons in the CC, Ncc, the EC, Nec, and in the ICR, 
NicR. Ncc = 4, Ncc = 2 with Nec = 2, Ncc = 3 
with Nec — Ij and Ncc ^ 2 with Nicr — 1. Since we 
do not use the muon system in eeee event reconstruction, 
we include events where the muon system was not fully 
operational. This leads to a slightly higher integrated 
luminosity in the eeee final state compared to the ee/i/x 
and fJ.fJ,fJ.fJ, final states. 



B. eefifj, channel 

The eeiiiJ, channel is divided into three sub-channels 
that depend on the number of electrons in the CC: 
Ncc = 2, Ncc = 1, and Ncc = 0. No ICR electrons are 
used in this channel. As in the eeee final state, we apply 
this splitting because the instrumental background con- 
tamination varies significantly depending on the number 
of central electrons. We require at least two electrons and 
two muons; if there are more leptons in the event, only 
the highest-pT leptons of each type are used. To reject 
cosmic ray background, the cosine of the angle between 
the muons must satisfy cosa < 0.96, and the acoplanarity 
[36] between the two muons must be greater than 0.05 
radians. We further require jAzuc^l < 3.0 cm between 
the muon tracks, where ZucA refers to the location along 
the beam axis where the track has its distance of closest 
approach to the beamline. Also, we impose the require- 
ment that Ai? > 0.2 between all possible electron- muon 
pairings. Both the muon pair and electron pair invari- 
ant masses must exceed 30 GeV. There is no opposite 
charge requirement placed on the lepton pairs in order 
to maximize acceptance. 



C. fifififi final state 

In the four-muon final state, there must be at least four 
muon candidates satisfying the requirements in Sec. IV, 
and at least two of the muons must be matched to tracks 
found in the muon system. The four-muon system must 
be charge neutral (Xli=i 9^ = ^)^ and only oppositely 
charged pairs are considered as Z boson candidates. If 
more than four muons are reconstructed in the event, 
we consider only the four highest-pT muons. We further 
require jAzcc"^! < 3.0 cm between all muons. We also 
require that one of the two possible sets of dimuons has 
a dimuon mass M^^ > 30 GeV for both dimuons. 



VI. INSTRUMENTAL BACKGROUND 

The instrumental background primarily arises from 
Z{-^ U) + jets and Z{-^ U)+'y + jets production (with 
smaller contributions from WZ + jets, WW + jets, W + 
jets, and multijct production with >4 jets). These events 
contaminate the four-lepton channels when a jet is falsely 
reconstructed as an isolated lepton. Z{-^ M) + -f + jets 
production where a photon and a jet are mis-identified 
as an electron contaminates the eeee and eefxfi channels. 

We estimate the instrumental background using the 
data. We first find the probability for a jet to be mis- 
identified as a lepton, Pj£. A tag and probe method is 
used to determine Pji where di-jet activity is considered 
with jet pt > 15 GeV. The tagged jet must be associated 
with a jet that fired a single jet trigger and be the highcst- 
Pt jet in the event. We then look for a probe jet with 
\A(j)\ > 3.0 with respect to the tag jet, where (j> is the 



azimuthal angle. To suppress contamination from VF-|-jet 
events, we require ]^rp < 20 GeV in the tag and probe 
sample. The probe jets form the denominator of the Pji 
calculation. 

To calculate the numerator of the Pje estimate, we 
first find all good electrons in the event with a p^ > 15 
GeV. We then select those electrons that satisfy the same 
criteria imposed on the probe jets, noted above. The Pje 
estimate is parametrized as a function of the jet pr and 

m- 

The Pjf^ estimate is determined using a similar 
method. The tagged jet is defined as was done for elec- 
tron events, but in the numerator, rather than have an 
electron, we use any muon that has |A0| > 3.0 from the 
tag jet, and Pj^ is taken as the number of muons di- 
vided by the number of probe jets in the sample. The 
Pj^ estimate is parameterized in terms of px and rj. 

The Pje estimates for both electrons and muons are on 
the order of 10^^. 

To estimate the instrumental background for the eeee 
final state, Pje is applied to events with three recon- 
structed electrons and one or more jet. The jet kinemat- 
ics are used to model the electron kinematics in the event. 
This method accounts for events where either a photon 
or a jet is misreconstructed as one electron and a jet 
is misreconstructed as the other. This method overesti- 
mates the background from events with two real electrons 
and two jets misreconstructed as electrons. To determine 
the rate, we look at events with two reconstructed elec- 
trons and two or more reconstructed jets and apply Pje 
to both jets. The number of ee plus two jet events after 
Pje is applied to both jets is found to be negligible, so 
only eee-fjet events are used to model the instrumental 
background distributions in the eeee final state. 

The instrumental background in the ee/x/x channel is 
calculated from two different contributions. The first 
contribution is from events with efifi plus one or more 
jet, where we apply Pje to the jet. This method gives 
an estimate of a background due to Z(— > iJ-ii) + jets and 
Z(— > iiii) + 7 + jets where a jet has been reconstructed 
as an electron. We also consider the ee plus two jet or 
more case, where we apply Pj^ to the jets. This method 
gives an estimate of the background due to Z{-^ ee) + 
jets where the jets can contain muons. 

The Pj^ is applied to jets in fi/i plus two or more jets 
data to determine the instrumental background for the 
fJ-IJ-fJ-fJ- channel. 

Background estimates derived from the above method 
can be found in Tables I-III in each final state. 



VII. SYSTEMATIC UNCERTAINTIES 

The following factors contribute to the systematic un- 
certainty on this measurement. We assess a 1% trigger 
efficiency uncertainty. Lepton identification uncertain- 
ties are calculated by studying Z ^ ££ events; lepton 
identification uncertainties of 3.7% per CC and EC elec- 



TABLE I: Contributions from non- negligible backgrounds in the eeee subchannels, plus expected t-channel ZZ and Higgs boson 
signals and number of observed events. Uncertainties are statistical followed by systematic. 





2 CC 
2 EC 


3CC 
1 EC 


4CC 


> 2 CC 
1 ICR 


Instrumental backg. 
Migration 


0.15 ±0.01 ±0.03 
0.014 ±0.001 ±0.002 


0.12 ±0.01 ±0.02 
0.023 ± 0.001 ± 0.004 


0.05 ±0.01 ±0.01 
0.025 ± 0.001 ± 0.004 


0.29 ± 0.04 to'.'ia 
0.024 ± 0.001 ± 0.003 


Total VLon-ZZ 
background 


0.17 ±0.01 ±0.03 


0.14 ±0.01 ±0.02 


0.08 ±0.01 ±0.01 


0.32 ± 0.04 to'^i 


Expected 

t-channel Z/7* Z/7* 


0.48 ±0.01 ±0.07 


1.14 ±0.01 ±0.17 


1.03 ±0.01 ±0.15 


1.47 ±0.01 ±0.19 


Expected gg ^f H 
Mh = 125 GeV 
Expected ZH 
Mh = 125 GeV 


< 0.001 
0.003 


0.001 
0.006 


0.004 
0.010 


0.002 
0.008 


Total Higgs boson 
Mh = 125 GeV 


0.003 


0.007 


0.014 


0.010 


Observed 
Events 





1 


2 


2 



TABLE II: Contributions from non-negligible backgrounds in the ee/i/i subchannels, plus expected signal and number of 
observed events. Uncertainties are statistical followed by systematic. 





CC 


1 CC 


2 CC 


Instrumental backg. 
tt 

Migration 
Cosmic rays 


0.11 ±0.01 ±0.03 
(0.2 If,'? ±0.6) X IQ-'^ 

(2.1 +«:? tU) X W-' 
< 0.001 


0.21 ±0.01 ±0.04 

(1.0 +o;;^±0.2) X 10^2 

(5.0 ± 0.8 +?;^) X 10--' 

< 0.003 


0.27 ±0.01 ±0.04 

(0.3 toA ± 0-3) X 10^2 

(4.8 tli ± 1.0) X 10--' 

< 0.006 


Total non-ZZ 
background 


0.12 ±0.01 ±0.03 


0.23 ±0.01 ±0.04 


0.27 ±0.01 ±0.04 


Expected 

t-channel Z/'y* Z/'j* 


0.43 ±0.01 ±0.06 


2.37 ±0.02 ±0.28 


4.13 ±0.03 ±0.49 


Expected gg -^ H 
Mh = 125 GeV 

Expected ZH 
Mh = 125 GeV 


< 0.001 
0.001 


0.002 
0.015 


0.007 
0.036 


Total Higgs boson 
Mh = 125 GeV 


0.002 


0.017 


0.043 


Observed Events 


2 


1 


2 



tron, 6% per ICR electron, and 3.2% per muon are used. 
There is a 10%-50% systematic uncertainty on the in- 
strumental background expectation in the various final 
states that is due to observed variations in Pj£ when 
changing selection requirements for the di-jet sample as 
well as limited statistics in the data samples used. We 
assign 20% uncertainty to the tt background. This cov- 
ers uncertainty on the theoretical production rate of 7% 
for mtop — 172 GeV [31], plus variation in the cross sec- 
tion due to uncertainty on the top quark mass, and also 



that on the rate at which the b quark from top quark 
decays is misidentified as an isolated lepton. We esti- 
mate a PDF uncertainty of 2.5% on all MC samples. We 
assign a 7.1% uncertainty on the ZZ cross section used 
to estimate the migration background and the ZZ back- 
ground to the Higgs boson search. A systematic uncer- 
tainty of 6.1% is assessed on the luminosity measurement 
[37]. We assess a systematic uncertainty on the ZZ pT 
distribution by reweighting the pythia ZZ pT to match 
a distribution derived from SHERPA MC [24]. The ZZ 



Pt systematic is between 1% and 7% for signal t-channel 
ZZ events, but has up to a 40% effect on the migration 
background. We also assess systematic uncertainties on 
the muon and electron energy resolution [38] , which lead 
to an uncertainty on the cross section measurements and 
Higgs boson production limits of less than 2%. For the 
Higgs boson search, we assess a theoretical uncertainty 
on the expected gluon fusion and ZH associated cross 
sections of 10.9% and 6.2%, respectively [27, 28]. 



TABLE III: Contributions from non- negligible backgrounds 
in the ^i^i^i^J, channel, plus expected i-channel ZZ and Higgs 
boson signal and number of observed events. Uncertainties 
are statistical followed by systematic. 





Number of Events 


Instrumental backg. 
Migration 
Cosmic rays 


0.12 ±0.01 to'.os 

(0.34 ±0.02 +'(5:oI) X 10"^ 

<0.01 


Total non- ZZ 
background 


0.15 ±0.011°;°^ 


Expected 

t-channel Z/'y* Z/'j* 


4.26 ± 0.02 ± 0.43 


Expected gg ^> H 
Mh = 125 GcV 

Expected ZH 
Mh = 125 GeV 


0.007 
0.033 


Total Higgs boson 
Mif = 125 GeV 


0.040 


Observed Events 


3 



VIII. CROSS SECTION MEASUREMENT 

The data are used to measure the production cross sec- 
tion pp — > ZZ at y/s = 1.96 TeV. The integrated lumi- 
nosities analyzed for the three channels are 9.8, 9.6, and 
9.6 fb~^ for the eeee, ee^^, and fJ-fJ-fifJ, channels, respec- 
tively. A summary of the signal and background event 
expectations are included in Tables I-III for the three 
channels. 

We observe five eeee candidate events, five eefxfi candi- 
date events, and three jjljijiji candidate events, for 13 data 
events total, with a total of 16.8 ± 1.9 (stat±syst+lumi) 
expected events. 

A negative log-likelihood function is constructed by 
taking as input the expected signal acceptance, the num- 
ber of expected background events, and the number of 
observed events in each of the subchannels. The signal 
acceptance times efficiency for each channel are shown 
in Tables IV-VI. The branching ratio for each chan- 
nel is calculated using the relevant Z boson branch- 
ing ratios from Ref. [39]. The cross section, a^ is 
varied to minimize the negative log-likelihood, which 
gives a{pp ^ Zh*Zh*) = 1.26l°1^ (stat)+°-;^ (syst) ± 
0.08 (lumi) pb for M(Z/7*) > 30 GeV. We then calcu- 
late the ratio of a{pp —)■ Z/"f*Z/j*) to a{pp -> ZZ) 
for this mass region using MCFM [25], and from this 
correction determine the pp — >■ ZZ cross section to be 
l-OStoJo (stat)+o;" (syst) ± 0.06 (lumi) pb. We combine 
this measurement with the pp — > ZZ cross section mea- 
sured in the i'^l~vv final state using data from the DO 
detector [11], giving a total combined pp — >■ ZZ cross sec- 
tion of 1.32l°J9 (stat) ± 0.12 (syst) ± 0.04 (lumi) pb. The 
measured cross section values are consistent with the SM 
expectation of 1.43 ± 0.10 pb [25]. 



IX. HIGGS BOSON PRODUCTION LIMITS 



The main Higgs boson production mechanisms that 
can result in four final state charged leptons are gluon 
fusion and ZH associated production. 

For Higgs boson events produced through gluon fu- 
sion, final states with four charged leptons arise from the 
decay H — )■ ZZ, where both Z bosons then decay lepton- 
ically. As all of the decay products of the Higgs boson 
in this decay are well measured, the best discriminating 
variable between the gluon fusion Higgs boson signal and 
the backgrounds is the four-lepton invariant mass. 

In the case of associated ZH production, two of the 
leptons in each event can come from the decay of the 
associated Z boson, so Higgs decay modes with two or 
more final state leptons will contribute to our signal. The 
majority of the ZH signal arises from H — >■ t~^t~ , H — > 
WW , and H — > ZZ decays. We expect large $rp in these 
events, due to the neutrinos from the t and W boson 
decays, as well as in events where one Z boson from the 
H — >• ZZ decays to neutrinos. 

We therefore set limits on SM Higgs boson production 



TABLE IV: Acceptance x efficiency for tlie eeee subcliannels, for ZZ — >■ eeee and ZZ — > eerr decays. Uncertainties are 
statistical followed by systematic. 



Channel 



2 CC, 
3CC, 
4CC 



2 EC 
1 EC 



eeee 



0.025 ± 0.001 ± 0.004 
0.059 ±0.001 ±0.011 
0.053 ± 0.001 ± 0.009 



> 2 CC, 1 ICR 0.076 ± 0.001 ± 0.012 0.0007 ± 0.0001 ± 0.0001 



0.0002 ±0.0001 ±0.0001 
0.0006 ± 0.0001 ± 0.0001 
0.0007 ±0.0001 ±0.0001 



TABLE V: Acceptance x efficiency for the ee/i/^ subchannels, for ZZ — >■ eejj,fi, ZZ — > eerr, and ZZ — > jijitt decays. 
Uncertainties are statistical followed by systematic. 



Channel 


ee^^ 


eCTT 


^flTT 


OCC 

1 CC 

2 CC 


0.011 ±0.001 ±0.001 
0.063 ±0.001 ±0.007 
0.110 ±0.001 ±0.012 


0.0001 ± 0.0001 ± 0.0001 
0.0007 ± 0.0001 ± 0.0001 
0.0014 ±0.0001 ±0.0002 


0.0002 ± 0.0001 ± 0.0001 
0.0007 ±0.0001 ±0.0001 
0.0019 ±0.0001 ±0.0002 



TABLE VI: Acceptance x efficiency for the pppp channel, 
statistical followed by systematic. 



for ZZ — > n^i^^Ji and ZZ — >■ fifirr decays. Uncertainties are 



lifififi 



0.224 ± 0.002 ± 0.022 



HflTT 



0.0032 ± 0.0002 ± 0.0003 



using the four-lepton invariant mass and the Ifrp. The 
four-lepton mass and ^rp are shown in Fig. 2, with the ex- 
pected Higgs boson signal distributions for a Higgs boson 
mass, Mh, of 125 GeV. Additional differential distribu- 
tions are provided in Appendix A. The expected yields 
for each production and decay mode for each Higgs bo- 
son mass considered are shown in Table VIL For events 
with $]j, < 30 GeV, the four-lepton mass is used to dis- 
criminate the Higgs boson signal from all backgrounds; 
in events with ]^rp > 30 GeV, the ^j- is used. For the 
Higgs boson search, the i-channel Z/j* Z/j* background 
is fixed to the SM expectation. 

We find no evidence of SM Higgs boson production 
and proceed to set limits. We consider potential Mh 
values between 115 and 200 GeV, in 5 GeV increments. 
We calculate limits on the SM Higgs boson production 
cross section using a modified frequcntist approach [40- 
42]. A log-likelihood ratio (LLR) test statistic is formed 
using the Poisson probabilities for estimated background 
yields, the expected signal acceptance, and the number of 
observed events for each considered Higgs boson mass hy- 
pothesis. The confidence levels are derived by integrating 
the LLR distribution in pseudo-experiments using both 
the signal-plus-background hypothesis (CLs+b) and the 
background-only hypothesis (CLi,). The excluded pro- 
duction cross section is taken to be the cross section for 
which the confidence level for signal, CLs = CLs+t/CLh, 
is less than or equal to 0.05. 

The calculated limits are listed in Table VIII. At 
Mh = 125 GeV, we expect to set a limit of 42.8 times 
the SM cross section at the 95% C.L., and observe a limit 



of 42.3 times the SM cross section. The limits vs. Mh 
are shown in Fig. 3, along with the associated LLR dis- 
tribution. 

X. CONCLUSIONS 



We have measured the production cross section for 
pp -^ Zji* Z/-f* with M(Z/7*) > 30 GeV to 
be 1.26+!]:^^ (stat);oJ5 (syst) ± 0.08 (lumi) pb. We 
correct this measurement by the expected ratio of 
(7{pp —>■ Z/7*Z/7*) to a{pp -^ ZZ) for this mass 
region and obtain a pp — > ZZ cross section of 
1.05l|]j^(stat)+o;j2(syst) ± 0.06 (lumi) pb. We also 
searched for the Higgs boson in the four lepton final state, 
assuming that the i-channel ZZ pair is produced with the 
cross section predicted by the SM. At Mh = 125 GeV, 
we expect a limit of 42.8 times the SM cross section, and 
set a limit of 42.3 times the SM cross section at the 95% 
C.L. 
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Appendix A: Differential Distributions 

Figs. 4-9 show differential distributions of the events 
used in the i-channel ZZ cross section measurement and 
Higgs boson search. Some of these distributions are kine- 
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FIG. 3: The (a) expected and observed 95% C.L. upper limits 
on the SM Higgs boson production cross section relative to 
the value expected in the SM, and the (b) log-likelihood ratio 
for all four lepton channels combined. 
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matic properties of dilepton systems; in the ee/x/i final 
state, the pairings of ee and fifi are always used. In 
the eeee and fJ.fJ,fJ,fJ, final states, there may be multiple 
combinations passing our selection requirements. If there 
are multiple passing combinations, we use the combina- 
tion that yields a dilepton pair with an invariant mass 
closest to the nominal Z boson mass of 91.2 GeV [39]. 
Fig. 4 shows the dilepton invariant mass and the px of 
the four-lepton system. The pt and r]d distributions for 
each lepton in our events are shown in Fig. 5 and 6, re- 



spectively. The Z/^* Pt distributions for the highest-p^ 
(leading) and second lepton pair are in Fig. 7. Fig. 8 
shows the distributions of the opening angles between the 
best matched lepton pairs in each event in the azimuthal 
angle, Aip, and AR. Fig. 9 shows the angle 0decay, which 
is the angle through which the lepton side of one of the 
Z/j* boson decay planes is rotated into the lepton side 
of the other Z/j* boson decay plane, and measured in 
the center-of-mass frame of the Z/Y* Z/j* system [43]. 
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FIG. 4: Distributions of (a) the dilepton invariant mass and (b) the transverse momentum of the four-lepton system in data, 
expected signal and background. There are two entries per event in the dilepton invariant mass distribution. The Higgs boson 
signal for Mr of 125 GeV is shown scaled by a factor of 40. 
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FIG. 5: Distributions of the transverse momentum in data, expected signal, and backgrounds for the (a) highest-pr, (b) second- 
highest-pT, (c) third- highest-pT, and (d) lowest-pT leptons in each event. The Higgs boson signal for Mh of 125 GeV is shown 
scaled by a factor of 40. 
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FIG. 7: Distributions of the Z/^* pt for the (a) leading and (b) second highest-pT lepton pairings in each event. In the eeee 
and lJ,fJ,fJ,fJ, channels, the combination shown is that with one dilepton mass most consistent with a Z mass of 91.2 GeV. The 
Higgs boson signal for Mh of 125 GeV is shown scaled by a factor of 40. 
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are two entries per event in both distributions. The Higgs boson signal for AIh of 125 GeV is shown scaled by a factor of 40. 
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FIG. 9: Shown in (a) is definition of (^dccay, adapted from Ref. [43] with permission. Shown in (b) is the distribution of the 
azimuthal 0dccay angle. In the eeee and jJ.fJ,jJ.fJ, channels, (^docay is calculated between the combination is most consistent with a 
Z mass of 91.2 GeV for one of the two dileptons. The Higgs boson signal for Mh of 125 GeV is shown scaled by a factor of 40. 



